Astronomy & Astrophysics manuscript no. sdBrot' paper single new 


©ESO 2012 


July 2, 2012 





Hot subdwarf stars in close-up view 
II. Rotational properties of single and wide binary subdwarf B stars * 

S. Geier 1 and U. Heber 1 



Dr. Karl Remeis-Observatory & ECAP, Astronomical Institute, Friedrich- Alexander University Erlangen-Nuremberg, Sternwartstr. 7, 
D 96049 Bamberg, Germany 



> 

l> 

ON 
NO 

NO 
O 
(N 



X 



Received Accepted 



(N 

o 

(N 

S3 
33 

I— ». 

on: 

(N 

6 

03 



ABSTRACT 



Subluminous B stars (sdBs) form the extremely hot end of the horizontal branch and are therefore related to the blue horizontal branch 
(BHB) stars. While the rotational properties of BHB stars have been investigated extensively, studies of sdB stars have concentrated 
on close binaries that are influenced by tidal interactions between their components. Here we present a study of 105 sdB stars, 
which are either single stars or in wide binaries where tidal effects become negligible. The projected rotational velocities have been 
determined by measuring the broadening of metal lines using high-resolution optical spectra. All stars in our sample are slow rotators 
(v rot sin i < lOkms -1 ). Furthermore, the v rot sin /-distributions of single sdBs are similar to those of hot subdwarf s in wide binaries 
with main-sequence companions as well as close binary systems with unseen companions and periods exceeding =i 1.2 d. We show 
that blue horizontal and extreme horizontal branch stars are also related in terms of surface rotation and angular momentum. Hot blue 
horizontal branch stars (r e g > 1 1 500 K) with diffusion-dominated atmospheres are slow rotators like the hot subdwarf stars located 
on the extreme horizontal branch, which lost more envelope and therefore angular momentum in the red-giant phase. The uniform 
rotation distributions of single and wide binary sdBs pose a challenge to our understanding of hot subdwarf formation. Especially the 
high fraction of helium white dwarf mergers predicted by theory seems to be inconsistent with the results presented here. 

Key words, binaries: spectroscopic - subdwarfs - stars: rotation 



1. Introduction 

Subluminous B stars (sdBs) show similar colours and spectral 
characteristics to main sequence stars of spectral type B, but are 
less luminous. Compared to main sequence B stars, the hydrogen 
Balmer lines in the spectra of sdBs are stronger while the helium 
lines are much weaker. The strong line broadening and the early 
confluence of the Balmer series is caused by the high surface 
gravities (log g 5.0 - 6.0) of these compact stars (/? s< ib — 0.1 — 
0.3 R Q ). Subluminous B stars are considered to be core helium- 
burning stars with very thin hydrogen envelopes and masses of 



about half a solar mass (Heber [T986l ) located at the extreme end 
of the horizontal branch (EHB). 

1.1. Hot subdwarf formation 
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The origin of EHB stars is still unknown (see Heber 2009 for a 
review). The key question is how all but a tiny fraction of the red- 
giant progenitor's hydrogen envelope was removed at the same 
time at which the helium core has attained the mass (^ 0.5 M Q ) 
to ignite the helium flash. The reason for this high mass loss at 
the tip of the red giant branch (RGB) is unclear. Several single- 
star scenarios are under discussion (D'Cruz et al. 1996; Sweigart 
[19971 De Marchi & Paresce [19961 Marietta et al. [20001 . which 
require either a fine-tuning of parameters or extreme environ- 
mental conditions that are unlikely to be met for the bulk of the 
observed subdwarfs in the field. 

According to Mengel et al. ( 1976), the required strong mass 
loss can occur in a close-binary system. The progenitor of the 
sdB star has to fill its Roche lobe near the tip of the red-giant 
branch (RGB) to lose a large part of its hydrogen envelope. 
The merger of close binary white dwarfs was investigated by 
Webbink (119841 1 and Iben & Tutukov (119841 1. who showed that 
an EHB star can form when two helium core white dwarfs 
(WDs) merge and the product is sufficiently massive to ignite 
helium. Politano et al. ([2008 ) proposed that the merger of a red 
giant and a low-mass main-sequence star during a common en- 
velope (CE) phase may lead to the formation of a rapidly rotating 
single hot subdwarf star. 

Maxted et al. (2001) determined a very high fraction of radial 
velocity variable sdB stars, indicating that about two thirds of 
the sdB stars in the field are in close binaries with periods of less 
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Fig. 1. T e g - log ^-diagram for the entire sample (not RV- 
variable) under study. The helium main sequence (HeMS) and 
the EHB band (limited by the zero-age EHB, ZAEHB, and the 
terminal-age EHB, TAEHB) are superimposed with EHB evo- 
lutionary tracks for solar metallicity taken from Dorman et al. 
(1993) labelled with their masses. Open circles mark objects 
where only upper limits could be derived for v rot sin i, filled cir- 
cles objects with significant v rot sin i. The size of the symbols 
scales with the value of v ro t sin i. 



than 30 days (see also Morales-Rueda et al. 2003 ; Napiwotzki et 
al. l2004al Copperwheat et al. [20TTb . Han et al. (120021 120031) used 
binary population synthesis models to study the stable Roche 
lobe overflow (RLOF) channel, the common envelope ejection 
channel, where the mass transfer to the companion is dynami- 
cally unstable, and the He-WD merger channel. 

The companions are mostly main sequence stars or white 
dwarfs. If the white dwarf companion is sufficiently massive, the 
merger of the binary system might exceed the Chandrasekhar 
mass and explode as a type la supernova. Indeed, Maxted et al. 
d2"000i ) found the sdB+WD binary KPD 1930+2752 to be a sys- 
tem that might qualify as a supernova la progenitor (see also 
Geier et al. 120071 ) . In Paper I of this series (Geier et al. [2010b) 
more candidate systems with massive compact companions, ei- 
ther massive white dwarfs or even neutron stars and black holes, 
have been found. Furthermore, Geier et al. (1201 ldl ) reported the 
discovery of an eclipsing sdB binary with a brown dwarf com- 
panion. 

1.2. Rotation on the horizontal branch 

The rotational properties of horizontal branch (HB) stars both 
in globular clusters and in the field all the way from the red to 
the blue end have been studied extensively in the last decades 
(Peterson [19831 [19851 Peters on et al. [19831 [19951 Behr et al 
I2000al |200 0bJ Kinman et al. [20001 Rec io-Blanco et al. 120021 
120041 Behr l2003all2003bl Carney et al. [20031 120081 . Most of 
these investigations were motivated by the puzzling horizontal 



branch morphologies in some globular clusters and the search 
for second or third parameters responsible for this phenomenon. 
The most interesting result of these studies is the discovery of a 
significant change in the rotational velocities of blue horizon- 
tal branch (BHB) stars when their effective temperatures ex- 
ceed ^ 11 500 K. HB stars cooler than this threshold value show 
v rot sin i values up to 40 km s , while the hotter stars rotate with 
velocities lower than lOkms" 1 . 

The transition in rotational velocity is accompanied by a 
jump towards brighter magnitudes in the colour-magnitude di- 
agram (Grundahl et al. 1999 ) and a change in the atmospheric 
abundance pattern. Stars cooler than ^ 11 500 K show the typical 
abundances of their parent population (e.g. For & Sneden l2010l ). 
while stars hotter than that are in general depleted in helium and 
strongly enriched in iron and other heavy elements such as ti- 
tanium or chromium. Lighter elements such as magnesium and 
silicon on the other hand have normal abundances (Behr et al. 
I2003all2003bl Fabbian et al. 120051 Pace et al. [2006t . Diffusion 
processes in the stellar atmosphere are most likely responsible 
for this effect. Michaud et al. ( 119831 ) predicted such abundance 
patterns before the anomalies were observed (see also Michaud 
et al. 120081 ). Caloi (Q999 ) explained the sharp transition between 
the two abundance patterns as the disappearance of subsurface 
convection layers at a critical temperature. Sweigart (2002]) in- 
deed found that thin convective layers below the surface driven 
by hydrogen ionization should exist and shift closer to the sur- 
face when the effective temperature increases. At about 12 000 K 
the convection zone reaches the surface and the outer layer of the 
star becomes fully radiative. Since convection is very efficient in 
mixing the envelope, diffusion processes do not operate in HB 
star atmospheres of less than 12 000 K. 

Slow rotation is considered as a prerequisite for diffusion. 
Michaud (1983 ) was the first to show that meridional circula- 
tion stops the diffusion process as soon as the rotational velocity 
reaches a critical value and could explain the chemical peculiar- 
ity of HgMn stars in this way. Quievy et al. (2009) performed 
similar calculations for BHB stars and showed that the critical 
rotational velocity is somewhere near 20kms~' at the transi- 
tion temperature of 1 1 500 K. This means that the atmospheric 
abundances of stars with lower v ro t sin i should be affected by 
diffusion processes. 

What causes the slow rotation that allows diffusion to hap- 
pen, is still unclear. Sills & Pinsonneault (2000 ) used a standard 
stellar evolution code and modelled the distribution of rotational 
velocities on the BHB. In order to reproduce the two popula- 
tions of fast and slow rotators they assumed two distinct main 
sequence progenitor populations with different rotational veloci- 
ties. In their picture the slowly rotating BHBs originate from 
slowly rotating main sequence stars. 

Another possible explanation is the spin-down of the surface 
layers by diffusion itself. Sweigart (2002) argued that the radia- 
tive levitation of iron triggers a weak stellar wind that carries 
away angular momentum. Vink & Cassisi (2002 ) showed that 
such winds are radiatively driven. 

Brown (120071 ) used a stellar evolution code including rota- 
tion and modelled the distribution of rotational velocities on the 
BHB. This code allows one to follow the evolution of the pro- 
genitor star through the He-flash. Brown (2007 ) argues that no 
significant angular momentum is exchanged between the stellar 
core and stellar envelope during the flash. The surface rotation 
of their models highly depends on the rotation of the surface 
convection zone, which contains most of the outer envelope's 
angular momentum. Hot BHB stars without surface convection 
zone rotate slower than the cooler ones with convection zone. 
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This approach allows one to reproduce the observed v 10t sin i- 
distribution of BHB stars without assuming bimodal stellar po- 
pulations (Brown et al. 2008). 

While the rotational properties of horizontal branch stars 
both in globular clusters and in the field are thoroughly exam- 
ined, the investigation of EHB stars has mostly been restricted to 
close binary systems, where tidal interaction plays a major role 
(Geier et al. 1201 0b). Very few apparently single EHB stars have 
been studied so far, all of which are slow rotators (< 10kms~', 
e.g. Heber et al. 120001 Edelmann I200TT ). 

In this paper we determine the projected rotational velocities 
of more than a hundred sdB stars by measuring the broadening 
of metal lines. In Paper I (Geier et al. [2010b) the rotational prop- 
erties of sdBs in close binary system were derived and used to 
clarify the nature of their unseen companions. Here we focus 
on the rotational properties of apparently single sdBs and wide 
binary systems, for which tidal interactions become negligible. 

In Sect. [2] we give an overview of the observations of high- 
resolution spectra and the atmospheric parameters of our sample. 
The determination of the rotational properties of 105 sdB stars 
are described in Sect. [3] the results are interpreted in Sect.|4]and 
compared to the corresponding results for BHB stars in Sect. [5] 
The implications for the sdB formation scenarios and the fur- 
ther evolution to the white dwarf cooling tracks are discussed in 
Sect. [6] and Sect. [71 respectively. Finally, a summary is given in 
Sect. [8] 

2. Observations and atmospheric parameters 

ESO- VLT/UVES spectra were obtained in the course of the ESO 
Supernovae la Progenitor Survey (SPY, Napiwotzki et al. 2001 , 
120011 at spectral resolution R =* 20 000 - 40 000 covering 3200 - 
6650 A with two small gaps at 4580 A and 5640 A. Ea ch of the 
50 stars was observed at least twice (Lisker et al. 2005). 

Another sample of 46 known bright subdwarfs was observed 
with the FEROS spectrograph (R = 48 000, 3750 - 9200 A) 
mounted at the ESO/MPG 2.2m telescope (Geier et al. l20T2t . 

Six stars were observed with the FOCES spectrograph (R = 
30000, 380 - 700 A) mounted at the CAHA 2.2m telescope 
(Geier et al. 120121 1. Two stars were observed with the HIRES 
instrument (R = 45 00 0, 360 - 5120 A) mounted at the Keck 
telescope (Heber et al. 2000). One star was observed with the 
HRS fiber spectrograph at the Hobby Eberly Telescope (R = 
30 000, 4260- 6290 A, Geier etal.EpJOB. 

Because a wide slit was used in the SPY survey and the 
seeing disk did not always fill the slit, the instrumental profile 
of some of the UVES spectra was seeing-dependent. This has 
to be accounted for to estimate the instrumental resolution (see 
Paper I). The resolution of the spectra taken with the fiber spec- 
trographs FEROS and FOCES was assumed to be constant. 

The single spectra of all programme stars were radial- 
velocity (RV) corrected and co-added in order to achieve higher 
signal-to-noise. 

Atmospheric parameters of the stars observed with UVES 
have been determined by Lisker et al. (12005b . HP 20 5805 and 
Feige 49 have been analysed by Przybilla et al. ( 120061 , the two 
sdB pulsators KPD 2109+4401 and PG 1219+534 by Heber et 
al. (120001 ). and the sdB b inaries PG 1 725+2 52 and TONS 135 
by Maxted et al. (120011 ) and Heber d 19861 ), respectively. The 
rest of the sample was analysed in Geier et al. (120121 ) and a 
more detailed publication of these results is in preparation. We 
adopted the atmospheric parameters given in Saffer et al. ( 1994) 
for [CW83] 1758+36. 



The whole sample under study is listed in TablesQ]and[2]and 
the effective temperatures are plotted versus the surface gravi- 
ties in Fig. [1] Comparing the positions of our sample stars to 
evolutionary tracks, we conclude that all stars are concentrated 
on or above the EHB, which is fully consistent with the theory. 
We point out that the inaccuracies in the atmospheric parame- 
ters do not significantly affect the derived projected rotational 
velocities. 

3. Projected rotational velocities from metal lines 

To derive v ro t sin/, we compared the observed spectra with 
rotationally broadened, synthetic line profiles using a semi- 
automatic analysis pipeline. The profiles were computed for the 
appropriate atmospheric parameters using the LINFOR program 
(developed by Holweger, Steffen and Steenbock at Kiel univer- 
sity, modified by Lemke[l997j). 

For a standard set of up to 187 unblended metal lines from 
24 different ions and with wavelengths ranging from 3700 to 
6000 A, a model grid with appropriate atmospheric parameters 
and different elemental abundances was automatically generated 
with LINFOR. The actual number of lines used as input for an 
individual star depends on the wavelength coverage. Owing to 
the insufficient quality of the spectra and the pollution with tel- 
luric features in the regions blueward of 3700 A and redward 
of 6000 A we excluded them from our analysis. A simultaneous 
fit of elemental abundance, projected rotational velocity and ra- 
dial velocity was then performed separately for each identified 
line using the FITSB2 routine (Napiwotzki et al. I2004bl ). A de- 
tailed investigation of statistical and systematic uncertainties of 
the techniques applied is presented in Paper I. Depending on the 
quality of the data and the number of metal lines used, an ac- 
curacy of about 1 .0 km s can be achieved. For the best spectra 
with highest resolution the detection limit is about 5.0 km s . 

Projected rotational velocities of 105 sdBs have been mea- 
sured (see Tables [T]|2]). Ninety-eight sdBs do not show any RV 
variability. In addition, seven are radial velocity variable systems 
with orbital periods of about a few days (see Table[2|). 

For eleven stars of our sample upper limits for the pro- 
jected rotational velocities have already been published (Heber 
et al. 2000; Edelmann et al. 2001) based on the same spectra 
as used here (see Table©. Only for PHL932 and PG 0909+276 
our measured v rot sin i deviate significantly from the results of 
Edelmann et al. (120011 ). This is most likely because they used 
fewer metal lines in their study. 

Przybilla et al. (2006) performed an NLTE analysis of 
Feige 49 and HD 205805 using the same FEROS spectra as we 
do here and derived a v rot sin i below the detection limit. Again 
our measurements are consistent with their results, because they 
are very close to the detection limit we derived for FEROS spec- 
tra of sdBs (=; 5 km s _1 , see Paper I). 

4. Projected rotational velocity distributions 

The projected rotational velocities of our full sample of 98 stars 
without radial velocity variations are all low (< lOkms -1 , see 
Table [I}. Taking into account the uncertainties, one can see that 
there is no obvious trend with the atmosperic parameters (see 
Fig.lD- 

Fig. [2] shows the distribution of v rot sin ; binned to the aver- 
age measurement error (1.5kms ). Eleven stars that had only 
fairly weak upper limits of 10 km s _1 , were sorted out. The 
distribution is very uniform and shows a prominent peak at 
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Fig. 2. Distribution of v rot sin i for the full sample. Objects with Fig. 4. Distribution of v 10t sin i for 16 sdBs with companions vis- 
limits below the detection limit have been stacked into the first ible in the spectra using the same binning as in Fig. [2] 
dotted bin. 




Fig. 3. Distribution of v rot sin i for 71 single stars from our sam- 
ple using the same binning as in Fig. [2] The solid grey line marks 
the distribution of v ro t sin i under the assumption of randomly 
oriented rotation axes and a constant v lot = 7.65 km s~', which 
matches the observed distribution very well. 



6-8 kms -1 . Because we can only determine the projected ro- 
tation, the true rotational velocities of most stars in the sample 
should beabout7-8 km s ~ 1 . 




0.5 



0.4 



§0.3 



0.2 



0.1 '<$•:■ 



10 15 
v rot sini [km s~'] 



20 



Fig. 5. Distribution of v ro t sin i for 8 radial velocity variable sdBs 
with orbital periods exceeding =* 1.2 d using the same binning as 
in Fig. |2] 



4.1. Single-lined sdBs 

Our sample contains 71 single-lined sdBs, of which the v rot sin i 
could be constrained. Ten stars of which we were only able to 
derive upper limits of 10km s were sorted out. Fig. [3] shows 
the v rot sin i distribution of this subsample. Most remarkably, 
the distribution is almost identical to that of the full sample. 
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Adopting a random distribution of inclination angles and a con- 
stant v ro t of =i 8kms _1 , the observed v rot sin /-distribution can 
indeed be well reproduced (see Fig. |2J. We therefore conclude 
that most single sdBs in our sample have very similar rotation 
velocities. 

4.2. Double-lined sdB binaries 

Our sample contains 18 sdBs with visible spectral signatures of 
cooler main sequence (MS) companions (e.g. Mgi, Lisker et al. 
2005 ). Again, two stars with upper limits of lOkms were ex- 
cluded. 

The orbital periods of these systems are long. Green et al. 
(2006) have argued that such systems should have periods of 
many months or years. Recently, Deca et al. (120121 ) were able 
to determine the orbital period P ^ 760 d of the sdB+K bi- 
nary PG 1018-047. Similar periods were reported by 0stensen 
& van Winckel ( 120121 ) for eight such binaries. The separations 
of the components are so wide that tidal interaction is negligible. 
Main-sequence companions do therefore not affect the rotational 
properties of the sdB stars in this type of binaries. 

The distribution for sdBs with composite spectra is displayed 
in Fig. [4] Taking into account the much smaller sample size, the 
result is again similar. We therefore conclude that the rotational 
properties of sdBs in wide binaries with MS companions are 
the same as those of single sdBs, although they have probably 
formed in a very different way (see Sect. |6]). 

4.3. Pulsating sdBs 

Two types of sdB pulsators are known. The slow pulsations of 
the V 1093 Her stars (sdBV s , Green et al. 120031 ) are not expected 
to influence the line broadening significantly (see Geier et al. 
2010b). For the short-period pulsators (V 361 Hya type, sdBV r , 
Charpinet et al. 119971 Kilkenny et al. 1997) unresolved pulsa- 
tions can severely affect or even dominate the broadening of 
the metal lines and therefore fake high v ro t sin i. Telting et al. 
(2008 ) showed that this happens in the case of the hybrid pul- 
sator Balloon 090100001 using the same method as in this work. 
Unresolved pulsations are also most likely responsible for the 
high line broadening (39kms _1 ) measured for the strong pul- 
sator PG 1605+072 (Heber et al. [19991 [2000b. 

Our sample contains three known long-period pulsators 
(PHL 44, Kilkenny et al. [20071 PHL 457, Blanchette et al. 
2008] LB 1516, Koen et al. l20T0b and tw o short-period 
ones (KPD 2109+ 4401, Billeres et al. [19981 PG 1219+534, 
O'Donoghue et al. [19991 . The v rot sin i of KPD 2109+4401 is 
indeed among the highest of all sdBs in our sample (10.5 + 
1.6km s~'), but it is unclear if this might not be partly due to 
unresolved pulsations. Jeffery & Pollacco (2000) measured RV 
variations of 2kms~' for KPD 2109+4401. Taking this into ac- 
count, the sdBs rotational velocity may be slightly lower than 
measured. The v rot sin / of the other pulsators are not peculiar. 

For most stars in our sample it is not clear whether they are 
pulsators or not, because no light curves of sufficient quality are 
available. Because only about 5% of all sdBs show pulsations 
detectable from the ground, one may conclude that the contami- 
nation by pulsators should be quite low. Thanks to the extensive 
photometric surveys for sdB pulsators conducted by Billeres et 
al. d2002b , Randall et al. d2006b and Ostensen et al. d20T0b , we 
know that 27 stars from our sample do not show short-period 
pulsations. 



Restricting ourselves to these objects and again excluding 
those with visible companions, we constructed a "pure" sample 
of 16 single sdBs, for which the rotational broadening is proven 
to be disturbed neither by the presence of a companion nor by 
pulsations. The associated v rot sin i distribution does not differ 
from the other distributions (see Figs. [2]|4](. We therefore con- 
clude that unresolved pulsations do not significantly affect our 
results. 

4.4. Radial velocity variable sdBs 

In Paper I we showed that the v rot sin i distribution of sdBs in 
close binary systems is strongly affected by the tidal interaction 
with their companions, but that this influence becomes negligible 
if the orbital periods of the binaries become longer than 1 .2 d. 
It is instructive to have a look at the v rot sin /-distribution of these 
long-period radial velocity variable systems. From Paper I we 
selected all seven binaries with periods longer than 1.2 d, for 
which tidal synchronisation is not established. We added the sys- 
tem LB 1516, a binary with yet unknown orbital parameters, but 
for which Edelmann et al. (2005 ) provided a lower limit for the 
period of the order of days[3. 

Fig. [5] shows the associated distribution. Given the small 
sample size and although two stars have somewhat higher 
v rot sin i = 10 - 12kms _1 , the distribution is again very simi- 
lar to the distributions shown before (see Figs.[2]|4]). Subdwarf B 
stars in close binaries obviously rotate in the same way as sin- 
gle stars or sdBs with visible companions if the orbital period is 
sufficiently long. 



5. Comparison with BHB stars 

Projected rotational velocities of BHB stars have been deter- 
mined for many globular cluster and field stars (Peterson et al. 
IT9951 Behr I2003al I2003bl Kinman et al. [20001 Recio-Blanco 
et al. 120041 ). The results are plotted against the effective tem- 
perature in Fig. [6] The characteristic jump in v rot sin / at a tem- 
perature of about ^ 11 500 K can be clearly seen. The sdB se- 
quence basically extends the BHB trend to higher temperatures. 
The v rot sin / values remain at the same level as observed in hot 
BHB stars. 

Comparing the v rot sin i of BHB and EHB stars, one has to 
take into account that the radii of both types of horizontal branch 
stars are quite different, which translates directly into very dif- 
ferent angular momenta. While sdBs have surface gravities log g 
between 5.0 and 6.0, the surface gravities of BHB stars range 
from logg = 3.0 to 4.0. The BHB stars with the same rota- 
tional velocities as EHB stars have higher angular momenta. 
Assuming rigid rotation, the same inclination angle of the ro- 
tation axis, and the same mass of ^ 0.5 M for BHB and EHB 
stars, one can calculate the quantity v rot sin ixg~ 1 ^ 2 , which is di- 
rectly proportional to the angular momentum. The surface grav- 
ities of the sdBs were taken from the literature (see Sect. [2]), 
those for the BHB stars from Behr (12003al [2003bt and Kinman 
et al. (12000b . Since Peterson et al. (119951 ) and Recio-Blanco et 
al. (|2C)04 ) did not determine surface gravities for their BHB sam- 
ple, we adopted a logg of 3.0 for stars with temperatures below 
^ 10 000 K and 3.5 for the hotter ones as suggested by the results 
of Behr (12003al I2003bb and Kinman et al. ©00'). 



1 TONS 135 was not included because the orbital period of ^ 4d 
given in Edelmann et al. ( 2005 1 is not very significant and shorter peri- 
ods cannot be excluded yet. 
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Fig. 6. Projected rotational velocity plotted against effective tem- 
perature. The grey squares mark BHB and some sdB stars taken 
from P eterson et al. ([19951 Beh r (12003a! I2003bl . Kinman et al. 
(2000), and Recio-Blanco et al. £2004). Upper limits are marked 
with grey triangles. The black diamonds mark the sdBs from 
our sample. The vertical line marks the jump temperature of 
11 500 K. 
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Fig. 7. v ro t sin i x g~ 1 ^ 2 plotted against effective temperature. The 
grey squares mark BHB and some sdB stars taken from Peterson 
et al. (TT995I Behr d2003al I2003bb . Kinman et al. (2000), and 
Recio-Blanco et al. (2004). Upper limits are marked with grey 
triangles. The black diamonds mark the sdBs from our sam- 
ple. The vertical line marks the jump temperature of 11 500 K. 
Typical uncertainties for the sdBs are given in the upper right 
corner. 



In Fig.|7]v 10t sin ixg~ 1 ^ 2 is plotted against T e ff. The transition 
between BHB and EHB stars is smooth. Since the progenitors of 
the EHB stars lost more envelope material on the RGB, the EHB 
stars are expected to have lower angular momenta than the BHB 
stars. This is consistent with what can be seen in the Fig. [7] 

6. Implications for hot subdwarf formation 

The uniform distribution of low projected rotational velocities 
in single and wide binary sdBs has consequences for the open 
question of hot subdwarf formation. As shown in this study, sdBs 
appear to rotate at low but spectroscopically detectable velocities 
of 8 - lOkms -1 . These results are remarkably similar to those 
derived for their cooler relatives, the BHB stars. Hot subdwarfs 
are likely formed through binary interaction or merging, which 
is also accompanied by a transfer of angular momentum. The 
rotational properties of sdB stars therefore allow one to constrain 
possible formation scenarios. 

6.1. Uniform rotation of EHB stars and mass loss on the 
RGB 

The rotational properties of sdBs residing on the EHB are very 
similar to those of hot BHB stars. The only exception is that the 
EHB stars obviously lost more envelope in the red-giant phase 
and therefore retained less angular momentum. How the enve- 
lope is lost does not affect the rotational velocities of sdB stars, 
since the v' rot sin /-distribution of RV variable systems with or- 
bital periods sufficiently long to neglect the tidal influence of the 
companion (Fig. [5]) is similar to those of apparently single sdB 



stars (Fig. [3]) and for sdB stars with visible main sequence com- 
panions (Fig. [4]). 

The abundance patterns of sdBs are dominated by diffusion 
processes very similar to those of the hot BHB stars (Geier et al. 
2010a). No surface convection zone should be present, and ac- 
cording to the model of Brown (120071 ) the angular momentum of 
the outer layers should be low. Stellar winds and magnetic fields 
may help to slow down the upper layers of the star. However, 
Unglaub (2008 ) showed that the weak winds predicted for sdB 
stars are most likely fractionated and are therefore not able to 
carry away the most abundant elements hydrogen and helium. 

Angular momentum gained or retained from the formation 
process may also be stored in the stellar core, which may be 
rapidly rotating. Kawaler & Hostler (2005) proposed such a sce- 
nario and suggested an asteroseismic approach to probe the ro- 
tation of the inner regions of sdBs. Van Grootel et al. (2008]) and 
Charpinet et al. (2008) performed such an analysis for the two 
short-period sdB pulsators Feige48 and PG 1336-018, respec- 
tively, and found no deviation from rigid rotation at least in the 
outer layers of these stars down to about half the stellar radius. 
But these results may not be representative, because both stars 
are in close binary systems and are synchronised by the tidal in- 
fluence of their companions (Geier et al. 1201 0b ). The rigid body 
rotation may have been caused by this effect and may not be a 
general feature of sdBs. Another setback of these analyses is the 
problem that p-mode pulsations are not suited to probe the in- 
nermost regions of sdBs. In contrast to that, g-mode pulsations 
reach the stellar core and it should be possible to measure the 
rotational properties of the whole stellar interior with asteroseis- 
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mic methods. With the availability of high-precision light curves 
from the Kepler and CoRoT missions, the analysis of g-mode 
pulsators became possible and first results have been published 
by van Grootel et al. ( 1201 Oi l and Charpinet et al. ( 1201 lbb . 

For the RV variable systems CE ejection is the only fea- 
sible formation channel. The systems with visible compan- 
ions may have lost their envelopes via stable RLOF. Very re- 
cently, 0stensen et al. (12012.1 ) and Deca et al. (120121 1 reported 
the discovery of sdB+MS binaries with orbital periods up to 
=* 1200 d, which may have been sufficiently close for mass trans- 
fer. However, the visible companions to the sdBs may still have 
been separated by too much for an interaction with the subd- 
warf progenitors. More detailed binary evolution calculations are 
needed to solve this problem. Common envelope ejection and 
stable RLOF form similar sdB stars, because in both cases the 
hydrogen envelope is removed and the helium burning should 
start under similar conditions. It would therefore not be surpris- 
ing if their v lot sin /-distributions were to look similar. 

6.2. Where are the He-WD merger products? 

The v rot sin /-distribution of the single sdB stars (Fig. [3} is par- 
ticularly hard to understand in the context of the WD merger 
scenario. If a certain fraction or even all of the apparently sin- 
gle sdBs would have been formed in this way, one would not 
expect a v rot sin /-distribution that resembles that of the post- 
CE or post-RLOF sdBs. Gourgouliatos & Jeffery (2006) showed 
that the merger product of two WDs would rotate faster than 
break-up velocity, if angular momentum were conserved. These 
authors concluded that angular momentum must be lost dur- 
ing the merger process. One way to lose angular momentum 
are stellar winds and magnetic fields. Another explanation may 
be the interaction with the accretion disc during the merger. If 
the less massive object is disrupted, it should form an accretion 
disc around the more massive component. The WD can only 
gain mass if angular momentum is transported outward in the 
disc. This process is expected to spin down the merger prod- 
uct (Gourgouliatos & Jeffery 2006 ). According to a model pro- 
posed by Podsiadlowski (priv. comm.), the merger is accompa- 
nied by a series of outbursts caused by the ignition of helium. 
These flashes remove angular momentum from the merged rem- 
nant and should slow it down to rotational velocities of less than 
20kms- 1 . 

However, even if it is possible to slow down the merged rem- 
nant of two He- WDs, it is very unlikely that the merger pro- 
ducts would have a v ro t sin /-distribution almost identical to sdBs, 
of which we know that they were formed via CE-ejection or 
maybe stable RLOF. This would require an extreme fine-tuning 
of parameters, unless there is an as yet unknown mechanism at 
work, which leads to uniform rotation of the radiative, diffusion- 
dominated atmospheres. It is therefore questionable whether our 
sample contains stars that were formed by an He-WD merger or 
a CE-merger event. If this is not the case and because of the size 
of our sample, it would be safe to conclude that the merger chan- 
nel does not contribute significantly to the observed population 
of single hydrogen-rich sdO/Bs in contrast to the models of Han 
etal. (120021 120031) . 

This conclusion is consistent with the most recent results by 
Fontaine et al. (120121 1. who studied the empirical mass distri- 
bution of sdB stars derived from eclipsing binary systems and 
asteroseismic analyses. The lack of sdB stars more massive than 
^ 0.5 M , which would be the outcome of the merger channel, 
led to the conclusion that mergers are less frequent in the forma- 
tion process of isolated sdB stars than predicted by theory. 



The only known single and fast rotating hot subdwarf star 
EC 22081-1916 (Geier et al. 120 11 al l may be the rare outco me of 
a CE merger event as suggested by Politano et al. (f2008 ). It is 
unique among ^100 sdBs of our sample. 

Possible candidates for WD-merger products are the helium 
rich sdOs (He-sdOs, Stroer at al. 120071 1. since Hirsch et al. 
(2009) measured v lot sin / values of 20 - 30kms~' for some of 
those stars. Although their velocities are not particularly high, 
they are significantly different from the typical v rot sin i of sdBs. 
However, while the He-sdOs were first considered as single stars 
(Napiwotzki et al. 2008 ), evidence grows that a fraction of them 
resides in close binaries (Green et al. 2008; Geier et al. 1201 let . 
At least those He-sdOs could not have been formed by a He-WD 
merger. 

6.3. Alternative formation scenarios 

Because the canonical binary scenario for sdB formation, which 
rests on the three pillars CE ejection, stable RLOF and He-WD 
merger, turned out to be very successful not only in explain- 
ing the properties of sdBs in the field (Han et al. 2002 2003), 
but also in globular clusters ( Han 2008) and the UV-upturn phe- 
nomenon in old galaxies (Han et al. 120071) . the possible lack of 
merger candidates poses a problem. 

Alternative formation scenarios such as CE ejection trig- 
gered by substellar companions (Soker 1998; Bear & Soker 
120121 ) may be responsible for the formation of apparently sin- 
gle sdBs. Evidence grows that such objects are quite common 
around sdB sta rs (e.g. Silvotti et al. 120071 Geier et al. 1201 ldl 
Charpinet et al. 1201 lab . In the light of the results presented here 
and other recent observational evidence, the conclusion has to be 
drawn that the question of sdB formation is still far from settled. 

7. Connection to white dwarfs 

Owing to their thin hydrogen envelopes, hot subdwarf stars 
will not evolve to the asymptotic giant branch (AGB-manque, 
Dorman et al. 1993). After about 100 Myr of core He-burning 
on the EHB and a shorter episode of He-shell burning, these ob- 
jects will join the WD cooling sequence. 

The rotational properties of single WDs are difficult to deter- 
mine. Owing to the high pressure in the dense WD atmospheres, 
the spectral lines of WDs are strongly broadened and hence do 
not appear to be suitable to measure v rot sin i. However, the H a 
line often displays a sharp line core, which is caused by NLTE 
effects. In a small fraction of the WD-population metal lines are 
visible. However, excellent high-resolution spectra are necessary 
to constrain the projected rotational velocity (Berger et al. 2005 ). 

The derived upper limits (^ 10 - 50kms~') are consistent 
with the much lower rotational velocities of pulsating WDs de- 
rived with asteroseismic methods (^ 0.2 - 3.5kms _1 , Kawaler 
2003). Most single WDs are therefore obviously rather slow ro- 
tators. The reason for this is most likely a significant loss of mass 
and angular momentum due to stellar winds and thermal pulses 
in the AGB-phase, as has been shown by Charpinet et al. (2009 ). 

The properties of WDs evolved from sdB progenitors on the 
other hand should be very different. Since the hot subdwarfs by- 
pass the AGB-phase, both their masses and their angular mo- 
menta are expected to remain more or less constant when evolv- 
ing to become WDs. 

The average mass of these sdB remnants (^ 0.47 M Q ) is ex- 
pected to be significantly lower than the average mass of nor- 
mal WDs (^ 0.6 M Q ). But more importantly, the rotational ve- 
locities of these WDs must be very high. We have shown that 
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single sdBs have small, but still detectable v rot sin i. Assuming 
rigid rotation and conservation of mass and angular momen- 
tum, the rotational velocity at the surface scales with the stel- 
lar radius. Because the radius decreases by a factor of about 
10, the rotational velocity should increase by a factor of about 
100. Assuming an average v lot 8kms _1 for single sdBs, WDs 
evolved through an EHB-phase should therefore have an aver- 
age v rot si 800kms _1 . Because about 1% of all WDs are ex- 
pected to have evolved through an EHB-phase, we expect a sim- 
ilar fraction of extremely fast rotating, low-mass WDs. These 
high v rot sin lvalues should be easily detectable even in medium- 
resolution spectra. The sample of WDs with observed spectra 
from the Sloan Digital Sky Survey (Eisenstein et al. 2006]) for 
example should contain more than 100 of these objects. 

8. Summary 

We extended a project to derive the rotational properties of sdB 
stars and determined the projected rotational velocities of 105 
sdB stars by measuring the broadening of metal lines using high- 
resolution spectra. All stars in our sample have low v rot sin i < 
lOkms -1 . For 75% of the sample we were able to deter- 
mine significant rotation. The distribution of projected rotational 
velocities is consistent with an average rotation of 8 km s" 1 
for the sample. Furthermore, the v rot sin /-distributions of sin- 
gle sdBs, hot subdwarfs with main sequence companions vi- 
sible in the spectra and close binary systems with periods ex- 
ceeding 1 .2 d are similar. The BHB and EHB stars are related 
in terms of surface rotation and angular momentum. Hot BHBs 
with diffusion-dominated atmospheres are slow rotators like the 
EHB stars, which lost more envelope and therefore angular mo- 
mentum on the RGB. The uniform rotation distributions of sin- 
gle and wide binary sdBs pose a challenge to our understanding 
of hot subdwarf formation. Especially the high fraction of He- 
WD mergers predicted by theory seems to be inconsistent with 
our results. We predict that the evolutionary channel of single 
sdB stars gives birth to a small population of rapidly rotating 
WDs with masses lower than average. 
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27 700 


13.0 B 


112 


0.71 


24 


7.2+1.0 


UVES 


KUV01542-0710 C 


27 800 


16.3 B 


58 


0.92 


8 


7.2 + 2.1 


UVES 


HE2156-3927 C 


28 000 


14.1 B 


62 


0.61 


16 


7.0 ± 1.2 


UVES 


EC 03591-3232 


28 000 


11. 2 V 


131 




34 


4.8 ± 1.0 


FEROS 


EC 12234-2607 


28 000 


13.8 B 


60 




19 


6.8 ± 1.4 


FEROS 


PG 2349+002 


28 000 


12.0 B 


68 




11 


5.7 ± 1.5 


FEROS 


HE2322-0617 c - no 


28100 


15.7 B 


62 


0.70 


15 


6.8 ± 1.3 


UVES 


PG0258+184 c - no 


28 100 


15.2 B 


48 


0.99 


12 


7.2+ 1.7 


UVES 


HE0136-2758" 


28 200 


16.2 B 


29 


1.20 


27 


< 5.0 


UVES 


HE0016+0044 no 


28 300 


13.1 B 


58 


0.67 


14 


6.5 ± 1.3 


UVES 


PG 1549-001 no 


28 300 


14.8 B 


45 


1.16 


20 


5.6 ± 1.1 


UVES 


HE 2349-3 135 


28 500 


15.6 B 


53 


1.13 


13 


10.0 ± 1.7 


UVES 


EC 01 120-5259 


28 900 


13.5 V 


73 




19 


5.8 ± 1.2 


FEROS 


HE0007-2212 no 


29 000 


14.8 B 


53 


0.64 


21 


7.4+1.0 


UVES 


LB 275* 


29 300 


14.9 B 


48 


1.16 


20 


5.6+1.1 


UVES 


EC 03263-6403 


29 300 


13.2 V 


32 




40 


< 5.0 


FEROS 


HE1254-1540 cno 


29 700 


15.2 B 


54 


0.75 


20 


7.2 ± 1.3 


UVES 


PG 1303+097 


29 800 


14.3 B 


51 




18 


6.1 ± 1.5 


FEROS 


HE 2222-3738 


30200 


14.2 B 


61 


0.83 


28 


8.7 ± 1.0 


UVES 


HE 2238- 1455 


30400 


16.0 B 


48 


0.80 


14 


< 5.0 


UVES 


EC 03470-5039 


30500 


13.6 V 


53 




9 


7.3 ± 2.0 


FEROS 


Feige 38 


30600 


12.8 B 


148 




34 


5.3 ± 1.0 


FEROS 


HE 1038-2326 c 


30600 


15.8 B 


34 


1.27 


28 


< 5.0 


UVES 


PG 1710+490 


30600 


12.1 B 


80 




11 


7.1 ± 1.6 


FOCES 


HE 0447-3654 


30700 


14.6 V 


44 




11 


7.3 ± 1.8 


FEROS 


EC 14248-2647 


31400 


12.0 V 


104 




14 


7.0 ± 1.5 


FEROS 


HE0207+0030" 


31400 


14.7 B 


27 


1.30 


7 


5.1 +2.3 


UVES 


KPD2109+440P 


31800 


13.2 B 


136 




9 


10.5 ± 1.6 


HIRES 


EC 02542-3019 


31900 


12.8 B 


65 




13 


7.3 ± 1.5 


FEROS 


[CW83] 1758+36 nb 


32000 


11. 1 B 


110 




5 


5.7 ± 1.4 


FOCES 


TONS 155 c 


32 300 


14.9 B 


35 


0.85 


14 


< 5.0 


UVES 


EC 21043-4017 


32400 


13.1 v 


65 




8 


5.6 ± 1.8 


FEROS 


EC 20229-3716 


32500 


11. 4 V 


153 




29 


4.5 ± 1.0 


FEROS 


HS2125+1105 C 


32500 


16.4 B 


29 


0.80 


8 


6.0 ± 2.4 


UVES 


HE1221-2618 C 


32600 


14.9 B 


35 


1.06 


11 


6.8 ± 1.6 


UVES 


HS 2033+0821"° 


32700 


14.4 B 


43 


1.14 


37 


< 5.0 


UVES 


HE0415-2417 no 


32 800 


16.2 B 


34 


0.83 


10 


< 10.0 


UVES 


EC 05479-5818 


33 000 


13.1 v 


81 




20 


5.8 ± 1.1 


FEROS 


HE 1 200-093 l cno 


33 400 


16.2 B 


30 


0.86 


12 


< 5.0 


UVES 
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System 


Teff 


m B 


S/N 


seeing 




v ro t sin i 


Instrument 




[K] 


[mag] 




[arcsec] 




[kms 1 ] 




PHL932 


33 600 


12.0 B 


102 


1.10 


12 


9.0 ± 1.3 


UVES 


HE1422-1851 C -™ 


33 900 


16.3 B 


14 


0.58 


10 


< 10.0 


UVES 


PHL555 


34100 


13. 8 B 


56 


0.88 


17 


6.9 ± 1.2 


UVES 


HE 1419-1205 c 


34200 


16.2 B 


28 


0.69 


16 


< 10.0 


UVES 


PG 1219+534 5 


34 300 


12.4 B 


140 




11 


5.7 + 1.4 


HIRES 


HS 2216+1833° 


34400 


13. 8 B 


54 


0.90 


11 


5.3 ± 1.6 


UVES 


HE 1 050-0630™ 


34500 


14.0 B 


59 


1.20 


28 


7.3 ± 1.4 


UVES 


HE 151 9-0708™ 


34500 


15.6 B 


20 


0.84 


8 


9.0 + 2.4 


UVES 


HE 1450-0957 


34600 


15.1 B 


32 


0.71 


6 


9.0 ± 2.4 


UVES 


EC 13047-3049 


34700 


12.8 V 


68 




5 


6.8 + 3.6 


FEROS 


HS 1710+1614™ 


34 800 


15.7 B 


38 


1.30 


13 


< 5.0 


UVES 


PHL 334 


34 800 


12.5 B 


87 




13 


< 5.0 


FEROS 


Feige 49 


35 000 


13.2 B 


119 




40 


6.2 + 1.0 


FEROS 


HE2151-100P 


35 000 


15.6 B 


42 


0.66 


6 


6.7 + 2.4 


UVES 


PG 0909+ 164 s 


35 300 


13.9 B 


52 




4 


< 10.0 


FEROS 


HE 1021-0255™ 


35 500 


15.3 B 


40 


1.61 


11 


< 10.0 


UVES 


PG 0909+276™ 


35 500 


13.9 B 


82 




13 


9.3 ± 1.4 


FOCES 


HE 0101-2707 


35 600 


15.0 B 


67 


0.85 


12 


8.1 ± 1.5 


UVES 


EC 03408-1315 


35 700 


13.6 V 


66 




11 


8.8 + 1.8 


FEROS 


HE1352-1827 C 


35 700 


16.2 B 


24 


0.85 


5 


8.2 ± 2.7 


UVES 


PG 1207-032™ 


35 700 


13.1 B 


50 


0.64 


9 


6.6 ± 1.6 


UVES 


HE 0019-5545 


35 700 


15. 8 B 


38 


0.76 


7 


5.9 + 2.3 


UVES 


GD619 


36 100 


13.9 B 


96 


0.81 


10 


6.1 ± 1.5 


UVES 


HE 1441-0558 c -™ 


36400 


14.4 B 


30 


0.70 


8 


6.9 ± 2.0 


UVES 


HE 0123-3330 


36600 


15.2 B 


48 


0.66 


8 


6.9 + 1.8 


UVES 


PG 1505+074 


37 100 


12.2 B 


153 




4 


< 5.0 


FEROS 


tttt 1 A i\h i Ammo 
Hh 140/+0033 


no (in 

37 300 


15.5 


35 


0.72 


9 


< 10.0 


UVES 


PG 1616+144™ 


37 300 


13.5 B 


44 




4 


< 10.0 


FEROS 


EC 00042-2737° 


37 500 


13.9 B 


37 




9 


< 10.0 


FEROS 


PHL 1548 


37400 


12.5 B 


90 




10 


9.1 + 1.6 


FEROS 


PB 5333™ 


40600 


12.5 B 


66 




2 


< 10.0 


FEROS 


[CW83] 0512-08 


38400 


11. 3 B 


124 




14 


7.7 ± 1.1 


FEROS 



Notes. The average seeing is only given if the spectra were obtained with a wide slit in the course of the SPY survey. In all other cases the seeing 
should not influence the measurements. c Main sequence companion visible in the spectrum (Lisker et al. 20051. s Pulsating subdwarf of V 361 Hya 
type. 'Pulsating subdwarf of V 1093 Her type. No short-period pulsations have been detected either by ™Billeres et al. (2002 ), '"Randall et al. 
( 2006 1 or ™0stensen et al. < f20T0l . 'Misidentified as CBS 275 in Lisker et al. J20051 . 



Table 2. Projected rotational velocities of radial velocity variable sdBs. 



System 


Teff 


m B /v 


S/N 


^lines 


V ro t sin ; 


Instrument 




[K] 


[mag] 






[km s-'] 




TONS 135 


25 000 


13.1 B 


45 


35 


6.4 ± 1.0 


FEROS 


LB 1516' 


25 200 


12.7 B 


58 


23 


6.0 ± 1.3 


FEROS 


PHL457 1 


26 500 


13.0 B 


59 


47 


6.1 ± 1.0 


FEROS 


EC 14338-1445 


27 700 


13.5 V 


71 


39 


8.9 ± 1.0 


FEROS 


PG 1725+252 


28 900 


11. 5 B 


45 


11 


7.4 ± 1.1 


HRS 


PG 1519+640 


30300 


12.1 B 


104 


11 


9.4+ 1.4 


FOCES 


PG2151+100 


32 700 


12.9 B 


69 


9 


9.0+ 1.7 


FEROS 



Notes. 'Pulsating subdwarf of V 1093 Her type. 
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Table 3. Comparison with literature. 



System 


This work 


Literature 


Reference 




v rot sin ( 


Vrot sin ( 






n -It 

[kms '] 


[kms '] 




KPD 2109+4401 


10.5 + 1.6 


< 10.0 


Heber 


PG 1219+534 


5.7 + 1.4 


< 10.0 


et al. 42000|l 


BD+48° 2721 


4.7 ± 1.4 


< 5.0 


Edelmann 


Feige 65 


7.2 ± 1.1 


< 5.0 


etal. J200l|l 


HD 205805 


4.5 + 1.0 


< 5.0 




HD4539 


3.9 + 1.0 


< 5.0 




LB 1516 


6.0+1.3 


< 5.0 




PG 0342+026 


6.2+1.0 


< 5.0 




PG 0909+276 


9.3 + 1.4 


< 5.0 




PHL 932 


9.0+1.3 


< 5.0 




Feige 49 


6.2+1.0 


0.0* 


Przybilla 


HD 205805 


4.5 + 1.0 


0.0* 


et al. (|2006|( 



Notes. 'Adopted value for line fits is below the detection limit. 
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